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^ ' Abstract. Subcluster interactions within clusters of galaxies produce shocks that accelerate non- 

thermal particles. We treat Fermi acceleration of nonthermal electrons and protons by injecting 
«vj . power-law distributions of particles during the merger event, subject to constraints on maximum 

particle energies. The broadband nonthermal spectrum emitted by accelerated electrons and protons 
is calculated during and following the subcluster interaction for a standard parameter set. The inten- 
sity of y-ray emission from primary and secondary processes is calculated and discussed in light of 

^ ' detection capabilities at radio and y-ray energies. 
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Rich clusters contain thousands of galaxies and are the largest gravitationally bound 
systems in nature. Masses for rich clusters are ~lO^^A/0, with ~5-10% of the mass 



Q^' found in a hot intergalactic gas at temperatures of 2-12 keV. Rich clusters emit thermal 



bremsstrahlung with luminosities L^ ~ lO"*"^ — 10"^^ ergs s^^[l]. Poor clusters contain 
^ I hundreds of galaxies, have total masses ^lO^^M©, have a hot intergalactic gas of tem- 

^ ■ peratures 1-5 keV and X-ray luminosities Lj^ ~ 10^^ — lO^'^ ergs s^' [2]. Approximately 

90% of the total mass of clusters is in the form of nonluminous dark matter. 

In the hierarchical merging cluster scenario, poor clusters merge together to form 
/\ • richer clusters. Approximately 30-40% of galaxy clusters show evidence of substruc- 

j^ ■ ture in both the optical [3] and X-ray wavelengths [4]. Velocity differences between the 

observed structures is ^1000-2000 km s^K With gravitational forces driving the inter- 
action between the two systems, cluster mergers are consistent with highly-parabolic 
orbits. Typical sound speeds within the intergalactic medium (IGM) are ^800 km s^\ 
so shocks will form at the interaction boundary of the two systems. Computer models of 
merging clusters support the development of shocks in the IGM [5, 6, 7]. Dimensional 
arguments show that a cluster merger releases 10^^ — 10^^ ergs of gravitational potential 
energy when initial separations are of the order ~Mpc. 

Only the most massive and X-ray luminous galaxy clusters have extended diffuse 
radio sources. With projected linear sizes ~1 Mpc, these diffuse sources have no known 
optical counterparts. The diffuse radio emissions have two distinct characteristics. The 
extended diffuse emission found in the central region of a galaxy cluster with a regular, 
azimuthally symmetric shape is known as a radio halo. The diffuse emission found on 
the cluster periphery are the cluster radio relics. These features often have irregular 



shapes with signs of filamentary structure. Radio relics are associated only with clusters 
that show evidence of a recent or ongoing merger event. 

Shock fronts that form in the IGM as a result of a cluster merger event are thought 
to be associated with the cluster radio relics. The shock compression will orient any 
existing cluster magnetic field into the plane of the shock. Radio relics are characterized 
by highly organized magnetic fields with field strengths in the ~1 /jG range with linearly 
polarized field lines in the vicinity of the shock [8]. The shock front will accelerate a 
fraction of the thermal particles within the IGM by first-order Fermi acceleration. 

Recent work has highlighted the importance of nonthermal radiation from particles 
accelerated by shocks formed in merging clusters. Loeb and Waxman [9] have proposed 
that cluster mergers are the dominant contributor to the diffuse y-ray background, pro- 
vided the efficiency to convert the available gravitational energy into nonthermal elec- 
tron energy is ~5%. Some unidentified EGRET sources are claimed to be associated 
with y-ray emission from galaxy clusters [10]. Excess EUV emission from Coma can be 
explained by nonthermal electrons accelerated at merger shocks [11]. Variations in radio 
surface brightness will result from superposition of cluster emissions [12]. 

In order to examine this question in more detail, we have modified a supernova 
remnant code [13] to calculate nonthermal radiation spectra from primary and secondary 
particles in merger shocks. We go beyond previous treatments by considering both 
primary electron and proton acceleration, a time-dependent treatment of radiation losses, 
and radiation signatures of secondaries from proton-nuclear interactions. 



MODEL 

We have adapted a supernova remnant code [13] to treat the cluster merger scenario. 
The code is designed to calculate time-dependent particle distribution functions evolving 
through adiabatic and radiative losses for electrons and protons accelerated by the first- 
order Fermi process at the cluster merger shock. 

The electron and proton distribution functions originate from a momentum power-law 
injection spectrum. In terms of total particle energy E = mjc^, the injection function is 
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where ;? = |3y is the dimensionless momentum and s is the injection index. The maximum 
particle energy Emax is determined by the maximum energy associated with the available 
time since the beginning of the merger event, by a comparison of the Larmor radius with 
the size scale of the system, and by a comparison of the energy-gain rate through first- 
order Fermi acceleration with the energy-loss rate due to adiabatic, synchrotron and 
Compton processes. Particle injection ceases after the age of the shock front exceeds 
tacc — 10^ yrs. The constant Q^ normalizes the injected particle spectrum over the 
entire volume V{t) swept out by the shock front, and is determined by 



Ei';'= / dt / dEEQ,,p{E,t)V{t). (2) 



We assume a total available energy E'f = Tje ^10^^ ergs, and an efficiency factor r\ep = 
5% for both protons and electrons. Although the injection index s depends upon the 
Mach number of the shock, here we present calculations for a fixed index s — 2. 

The time evolving particle spectrum is determined by solving the Fokker-Planck 
equation in energy space for a spatially homogeneous IGM, given by 

The quantity Etot{E^t) represents the total synchrotron, Compton, Coulomb, and adia- 
batic energy-loss rate for electrons, and the sum of the Coulomb and adiabatic energy- 
loss rates for protons. Both protons and electrons are subject to diffusion in energy 
space by Coulomb interactions. The protons experience catastrophic losses due to 
proton-proton collisions on the time scale ipion- The spectra of secondary electrons and 
positrons are calculated from pion-decay products, and are subject to the same physical 
processes as the primary electrons. 

The synchrotron, Compton, bremsstrahlung, and pion-decay y-ray spectral compo- 
nents are calculated from the particle spectra following the methods described by 
Stumer, et al. [13]. We use a standard parameter set with a mean IGM number density 
niGM = 10^ cm^, a uniform cluster magnetic field of 5 = 0.1 ^uG, a constant shock 
speed Vs- = 1000 km s^\ and an acceleration period of 1 Gyr. Thus, V{t) = A{t)vst. We 
assume a constant surface area with a 1 Mpc radius. Particles loose energy due to adia- 
batic expansion according to the relation —y/j— V{t)/V{t) = l/t. 



RESULTS AND DISCUSSION 

Figure 1 shows nonthermal photon spectra calculated at 1 Gyr and 5 Gyr for our stan- 
dard parameters. The system is very luminous at radio frequencies during the particle 
acceleration phase because of intense Compton losses of nonthermal electrons on the 
CMB. After the acceleration period ends, the radiation from primary electrons is dom- 
inated by emission from secondaries. The n^ bump at 70 MeV is hidden beneath the 
Compton- scattered CMB radiation from primary electrons during the acceleration pe- 
riod, but dominates the y-ray spectrum after the acceleration period is over. The inten- 
sity of the n^ bump depends sensitively upon the relative efficiencies ripg for proton and 
electron acceleration. If no pion signature is found in nonthermal y rays from merging 
clusters, then Tj^ ^ Tje for our standard parameters. 

From the total photon spectra in Figure 1, we have calculated light curves at vari- 
ous photon energies in Figure 2 for a merging cluster at a distance of 100 Mpc. The 
maximum peak 1.4 GHz radio flux density of ~ 12 Jy occurs at the end of the acceler- 
ation period. The radio emission from the merger event is, however, distributed over an 
angular region ~0.5° of the merger shock, making it difficult to detect. 

Figure 2 presents calculations of the energy fluxes from a merger shock at 100 Mpc, 
for photon energies > 100 MeV, >1 GeV, and >300GeV. The maximum luminosity 
occurs at the end of the acceleration period due to the accumulation of nonthermal 
protons. Cooling times due to synchrotron and Compton energy losses are very short for 
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FIGURE 1. Total nonthermal photon energy spectra from a cluster merger shock. The solid curves are 
the total photon spectra summed from all processes. The dashed curves are the photons produced from 
the primary electrons, and the dotted curves are the secondary emission components. The left and right 
panels show the total nonthermal photon spectra at f = 1 Gyr and / = 5 Gyr, respectively, after the onset 
of the cluster merger event. 



high-energy electrons; electrons that Compton- scatter microwave photons to energies 
E = lOO^ioo MeV will cool on a time scale of tcooi — 6 x lO^/^/E[oo yrs. Consequently, 
the y-ray emission declines sharply after particle acceleration ceases, and approaches a 
level where all the y-ray production originates from proton interactions. The maximum 
flux for secondary electrons is ^2 orders of magnitude less than the maximum flux for 
the primary electrons, as can be seen by comparing the catastrophic proton loss time 
scale of ~30 Gyr with the 1 Gyr injection time scale. Because of adiabatic losses and 
the long catastrophic pion production cooling times, protons inject a slowly declining 
rate of secondary electrons in the 1-4 Gyrs after the initial acceleration event. 



SUMMARY 



We have studied the acceleration of protons and electrons in shocks formed by merging 
clusters of galaxies, and calculated the expected synchrotron, Compton, bremsstrahlung, 
and pion emission from the accelerated particles. Both the radio flux density at 1.4 GHz 
and the y-ray flux at photon energies > 100 MeV, > 1 GeV and >300 GeV for a source 
at 100 Mpc as a function of time are presented. The flux is greatest at the end of the 
particle acceleration phase, and merging clusters of galaxies at a characteristic distance 
of 100 Mpc should be detectable with radio telescopes and GLAST. Radio and y-ray 
detectability of merger events must contend with the angular extent of the emission. 

Because the merger event lasts for only ~ 1 Gyr, a large fraction of clusters will no 
longer be experiencing ongoing nonthermal particle acceleration from cluster mergers. 
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FIGURE 2. Calculated 1 .4 GHz radio flux density (left panel) and integrated flux of high energy 
photons (right panel) from a cluster merger shock. Both panels are normalized for a system at a distance 
of 100 Mpc. The radio flux density is in the 1.4 GHz band, and the integrated energy fluxes for the high 
energy photons are shown in the right panel. Also shown is the >300 GeV emission for a merger where 
Vj — 2000 km s^' and B = I /jG as the boldfaced, long-dashed line. 



GLAST, with a limiting sensitivity of 4 x 10^^-^ ergs cm^^ s^^ for a one year all- 
sky survey, will still be able to detect such systems if they are at optimal distances. 
Depending on the shock speed and magnetic field, cluster merger events can also be 
detectable sources of TeV radiation (see Fig. 2). We therefore expect detectable y-ray 
emission from clusters at distances ~100 Mpc. Within this distance a number of radio 
relics have been detected [14, 15]. 
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